Phantom scalar fields result in inflation rather than 

Big Rip 



There exists a variety of exact solutions of the scalar field Einstein equations, 
allowing for "phantom regions" with the negative kinetic field term. These regions 
can be cut out, their boundaries being sewn together in such a way that neither the 
scale factor (along with its first two derivatives) nor density or pressure will experi- 
ence a jump. Such a domain surgery eliminates the "Big Rip" scenario, substituting 
for it the standard inflation. 



1 Introduction 

The paper [1] suggested a new concept of death of the universe, called the "Big Rip" . Its 
main idea is to take into account the new hypothetical dark energy components with the 
state equation p = wp, with w < — 1. The matter with this equation-of-state parameter 
is called "phantom energy" [2] 

For the flat FLRW universe this yields the scale factor as a function of cosmic time 
(in the synchronous reference frame) as follows 1 : 



where a 2 = — 1 — w, to is a fixed time instant, and p(t) is the density. Clearly as t — > 
the scale factor blows up, underlying the Big Rip. 

The scalar field, minimally connected with gravitation and describing the hypothetical 
phantom component of dark energy has a negative kinetic term. Quantum theory of 
such fields is studied in [3], see also the references therein. Analysis has shown that 
fundamental phantom fields are an impossibility. However effective models of this sort, yet 
problematic, should not be ruled out. For example, such terms may arise in supergravity 
[4], in higher-derivative-gravity theories [5]; see also stringy phantom energy [6] and brane- 
world phantom energy [7]. 

In this note we attempt to demonstrate that the Big Rip can apparently be ruled out 
for purely classical (rather than quantum) reasons. There exists a class of exact solutions 

Everywhere in the sequel the common in cosmology system of units with 8irG = c = 1 is used. 
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of the Einstein equations, allowing the "phantom regions" . It will be further shown that 
all these regions can be cut out, with the edges being newly sewn up in such a way that 
neither the scale factor (and its first two derivatives) nor density or pressure experience 
a jump. 

2 Main idea 

Consider the simplest flat universe (k = 0) model with the FLRW metric: 

ds 2 = dt 2 - a(t) 2 (dx 2 + dy 2 + dz 2 } , 
filled with a scalar field with the Lagrangian 

L = ±crd lt <j>d v <f>-v(<f>). 

It's easy to verify that in the case when = 4>(t), the Einstein equations with the 
cosmological term 

are reduced to the Schrodinger equation 

^ = 3(V + A)i;, (1) 

with ip = a 3 . It is assumed that V = V(<fi(t)) = V(t), while the field time dependence is 
described by the following equation: 

Super-inflation (i.e. the case of dominating phantom fields) comes into play if H > 0, i.e. 
when a > a 2 /a. Suppose the function tp(t) has a point singularity at i*: 

Then it follows from (2) that in some neighborhood of t* the kinetic term becomes neg- 
ative, i.e. there exists an instant t — t\ < t* when 0(ti) = 0. Then t > t\ marks the 
phantom region, where the kinetic term is negative. To this effect, there are two possibil- 
ities: either this region stretches up to infinity into the future, or there exists an instant 
t 2 to the right of the singularity t*, such that 

<P(h) = 4>{t 2 ) = 0, 4> 2 {t) < 0, h < U < t 2 . 

It is easy to convince oneself that the latter alternative takes place if (i) the cosmological 
constant is zero and (ii) the self-action potential fairly quickly (for instance, as some 
power, see the next section) vanishes as t — > oo 2 . Indeed, H Solves the Ricatti equation 

H = V(t) + A - 3H 2 , 

2 More generally one may have V(t) + A — > fast enough as t — > oo, but such a condition implies 
adjusting the cosmological constant A = — V as , which hardly appears to be reasonable. Presumably, it 
is more sensible to render A = 0, the more so as that as of today, there would be no explanation for a 
small, but nonzero value of A. 



that is for A = 0, V — > one has H — > —3H 2 < 0, i.e. the phantom energy domination 
region is not present over large times. Then the quantity <fi must possess another zero, 
corresponding to some t 2 > q.e.d. 3 

The above conditions being satisfied, the universe would generally speaking enter the 
inflation phase as t — > t\. Indeed, let us assume that it is born at some time instant 
to < t\ and suppose there exists a time interval, during which all the energy conditions 
are satisfied. One can demonstrate that as time goes to t±, it is the strong energy condition 
to be violated first, the weak one being extant. As a result, the universe would get inflated 
up to the instant t±. Furthermore, when t — t± we arrive precisely at De Sitter's phase, 
but only instantaneously. Further follows the "phantom region", which should be cut 
out by matching the solutions corresponding to t — t\ and t = t 2 . Note that gluing 
the solutions in such a way should take place under constraints more stringent than 
those characteristic of similar problems of non-relativistic quantum theory based on the 
Schrodinger equation. The discussed problem calls for a C 2 gluing, namely the functions 
are to be matched up to their second derivatives. The latter requirement is to eliminate 
the possibility of discontinuities of pressure and density. A general solution of equation (1) 
is characterized by two integration constants, and it is not obvious that three matching 
conditions can be satisfied. Nevertheless, in is easy to verify that the third condition, 
namely the equality of the second derivatives, is automatically fulfilled! provided that 
the first two, namely the C 1 matching has been done. This apparent fluke is due to the 
fact that the matching is done at an inflection point of log^. 

The outcome is a solution without jumps of the scale factor and its first two deriva- 
tives, the density and its first derivative, and the pressure. Higher order derivatives of a 
(starting form the third one) p (starting from the second) and p (starting from the first 
derivative) do have jumps however. Pressure and density are connected via an algebraic 
state equation, and therefore the evolution of p(t) can be described in terms of that of 
p{t). In turn, p{t) = 3a{t) 2 / a{t) 2 . Thus the whole evolution can be computed from the 
state equation on the basis of the equation of motion 

— ^ + 3p(a, ct)^ a. 

This equation is strongly non-linear, but is nevertheless the second order ODE determined 
by the initial data. As one can see, the discontinuities pertain to higher order derivatives, 
which do not enter the above equation. Thus this equation plus the state equation do 
determine the dynamical evolution of the universe unambiguously, without the region of 
phantom field domination. 

Therefore, the existence of phantom regions entails cosmologies with standard infla- 
tion, rather than the universal Doomsday. If all along the solutions possess the right 
asymptotics (e.g. a Two thirds law dynamics) then the universe would automatically exit 
the inflationary phase without any extra arrangements. Let us further illustrate it by an 
easy example. 




3 Note that the exit from the phantom region would also take place for negative values of the cosmo- 
logical constant and a vanishing potential. No fine tuning is required in this case. However (mayhap 
erroneously!) we prefer to render A = 0. 



3 Model 



The law a ~ t 2//3 results in a potential 



The general solution of equation (1) with potential (3) and a zero cosmological constant 
is 

^ = ~j- c 2 t\ (4) 

where > are integration constants 4 . The solution describes two different universes, 
see Fig. I 5 . The first universe stretches into an unlimited past and has a finite order 
singularity at a point t = — {ci/ci) 1 ^. The second one looks rather mysterious: it 
emerges from a singularity t = to, then expands and experiences super-inflation, that 
is by the time t — the scale factor blows up, causing the Big Rip. Further follows a 
contraction, with a transition to a finite expansion, the scale factor changing the sign. 
Naturally this solution cannot be taken immediately for what it is. That is the super- 
inflation phantom region (to the left and to the right of the instant t — 0) should not be 
allowed, as the weak energy condition is violated therein. It is not difficult to verify that 
the pressure and density have the form 

ci (ci + 4c 2 t 3 ) (2c 2 t 3 - Cl f 

P ~ t 2 ( Cl + c 2 t 3 ) 2 ' P ~ W{ Cl + c 2 tZ) 2) 

and therefore the phantom region is confined to an interval from ti to t 2 , where 
Clearly, the solution on the right has to be determined by other integration constants: 

V (+) = y + M 2 . 

One can show that in order to sew up the solution consistently, one should take 

6i = (3 - 2V2) ci > 0, b 2 = (3 + 2V2) c 2 > 0. 

However, it is not entirely satisfactory as ti 7^ t 2 . It is not difficult to do a time shift 
within tp^, in order to ensure the following two conditions: 

1. The initial singularity corresponds to t — 0. 

2. ti = t 2 , i.e. there is no time discontinuity. 

Skipping the calculation, here is the answer: 
1) For < t < ti the dynamics is described as follows: 

= -TTT7 - 4 {t - Cl/C2)2 • l ' M = T, ^ r 

t - C1/C2 (t - C1/C2) 



4 The sign of V ( an d scale factor) is inessential: if a < then one must multiply it by —1. This 
because ip(t) is a solution of linear equation. If t = to is the zero of function i[)(t) and ip(t < t ) < and 
ij){t > t a ) > (or vice versa) then we have singularity at t = to so ip(t < 0) and i[)(t > 0) describes two 
different universes. 

5 All the figures are relegated to the Appendix. 



2) For ti < t it is 

^ {+ \t) = C ^ 3 + 2V ~ 2) + 4(3 - 2V2) (t - gf , V™ = 

t-g it- 9) 

where 



2 • 



For convenience, a pair of new constants c\ 2 has been introduced; also t 2 = t\. Matching 
has already been done: 



pW(^), p<>-Hti) = p (+) (ti) = o, = p (+) (ti). 

The graph a(t) for the second universe is given in Fig. 2. Note that for t < t\ the 
curve is described by the formula = (tp^ (t)) 1 ^ 3 , and for t > ti. = (ip^^t)) 1 ^ 3 
( c i — C2 = 1). Fig. 3 shows the time dependence of the acceleration a. The kink is 
located at the instant t = t\. 

As one can see, inflation over a finite interval takes place (in the region jutting over 
the t axis): if t < t\ then 

/c 7 \ 1/3 

ci(f -> ti) -> 0.127 Ml >0. 
Note that as t — > 00 the potential behaves as 

V~\e~^+, 0~-^=logt. 

Therefore, having removed the phantom region, we have arrived at a simple inflation 
model with an exit and asymptotic transition to Friedmann's universe expansion regime 
in the phase of the domination of matter. Clearly, the example above does not describe 
our Universe, as it lacks the oscillation phase. Nor it has the secondary acceleration. 
Nevertheless, it clearly indicates that one can use the suggested approach in order to 
produce in principle model cosmologies with exit, avoiding the parameter adjustment. 
Nor the hybrid models appear to be required. 

Let us recall that the solution describes one more, contracting universe, without a 
beginning but with an end to it. Solutions of this type look rather exotic. However, 
one inevitably runs into them, as the equations of General relativity allow changing the 
direction of time. 



4 Conclusion 

It has been shown above that the appearance of phantom regions does not entail super- 
inflation and the Big Rip (at least o the classical level). On the other hand, the existence 
of these regions is quite useful for constructing inflationary cosmologies with exit, without 
parameter adjustments. 



It does not follow that super-inflation is absolutely impossible. However, a naive 
approach to phantom regions does not appear to suggest the correct scenario for evolution 
of the Universe. As has been shown, a minor "surgery" eliminates the phantom regions, 
at least in the context of self-action potentials vanishing at infinity, in the case of a 
non-positive cosmological constant. 
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